Green's function approach to the magnetic properties of the kagome antiferromagnet 
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The S — 1/2 Heisenberg antiferromagnet is studied on the kagome lattice by using a Green's 
function method based on an appropriate decoupling of the equations of motion. Thermodynamic 
properties as well as spin-spin correlation functions are obtained and characterize this system as 
a two-dimensional quantum spin liquid. Spin-spin correlation functions decay exponentially with 
distance down to low temperature and the calculated missing entropy at T = is found to be 
0.46 In 2. Within the present scheme, the specific heat exhibits a single peak structure and a T 2 
dependence at low temperature. 

PACS numbers: 75.10. Jm,75.40.-s, 75. 50.Ee 



I. INTRODUCTION 

Antifcrromagnetic spin systems on fully frustrated lat- 
tices show mara/|-unusual behaviors in magnetic and ther- 
mal properties™! One of the main ingredients is that 
their unit cell allows for a continuous degree of free- 
dom and that the connectivity (corner sharing) allows 
for an extensive number of these degrees of freedom in 
the classical ground state. More subtle phenomena ap- 
pear when looking at quantum models on these lattices as 
quantum fluctuations may induce very unsusual ground 
states. In particular, wether the quantum ground state 
on fully frustrated lattices may not break the lattice sym- 
metry nor the spin group symmetry is still a highly de- 
bated question. A lot of work has been done to an- 
swer this question in the framework of the proposition 
of Andersond, looking first at Resonating Valence Bond 
states. 

One of the most studied candidates is the quantum 
5=1/2 Heisenberg antiferroma | gjoet- | pn | the , kagome lat- 
tice. Using various methodsQBltO'QBliaoll3, it has been 
shown that the low temperature physics should be dom- 
inated by short range RVB states which produce a con- 
tinuum of singlet states between the 5 = ground state 
and the first excited S = 1 state. A still controversial 
question is the presence of a very low temperature peak 
in the specific heat, much below the one corresponding 
to the onset of short range correlations, that could be 
ascribed to a high density of singlet states in the singlet- 
triplet spin gap. Associated to this low temperature peak 
is the missing (or not) entropy that would characterize 
an ordered or a disordered ground state. 

The experimental relevance of the model is extremely 
fragile as, in real compounds, many other parame- 
ters may drive the physics to very different universal- 
ity classes, described thijaugh inclusion of next near- 
est neighbcp-i interactions!!^, disordeni-J, antisymmetric 
interactions£3, etc... Nevertheless, in many cases, the 
Heisenberg antiferromagnet is a good starting point as 
it is generally believed that many properties of realistic 



systems come from deviations from the HeisenbereJi 
Examples are the layered oxide SrCrgpGa^-gpOigEjII 
(5 = 3/2) or._the organic pseudo-5 = 1 compound m- 
MPYNN.BF 4 Ei 

In this paper, the properties of the Heisenberg 5=1/2 
kagome antiferromagnet are addressed by using a spin 
Green's functions technique. One of the advantages of 
the method is that it is well suited for magnetic systems 
with no long range order as it uses a decoupling scheme 
based on short ranged spin correlations. This method 
was previously introduced by Kondo and Yamag *3 in 
the context of the one-dimensional Heisenberg model 
as a theory of spin-waves in absence of long-range or- 
der. More_xecently, it has been used to address different 
problepsBB and can be extended to include magnetic 
phasesEl Here, the formalism is used to compute ther- 
modynamic and magnetic quantities at all temperatures, 
like the internal energy, specific heat, entropy, magnetic 
susceptibility and structure factor. Correlation functions 
are also computed at all temperatures and for various 
separation distances. In the next section, the approxi- 



mation is presented and results are discussed in Sec. Ill 



II. MODEL AND APPROXIMATION 

The antiferromagnetic Heisenberg model is defined by 
the Hamiltonian 



H — Jij Si ■ Sj 



(1) 



where we assume 5 = 1/2 and a nearest-neighbor ex- 
change 



J(> 0) if i,j are nearest neighbors 



1 otherwise. 
The spin susceptibility 

Xij (u) =<5r;5+> Ll 



(2) 



(3) 
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FIG. 1: The kagome lattice. 

is obtained by a Fourier transformed Green's function 
>o,= J e lwt dt (4) 

where 

« A; B » t = -iO(t) ( B(0)] ) . (5) 

^From the equation of motion of the operator A(t) 
in the Heisenberg representation, the Green's functions 
above must satisfy the equation 



u < A; B > w = ([A, B]) + < [A, ft] ; B > a 



(6) 



The spin susceptibility of the present model can be 
obtained by 



wxjj-(w) = -2% (Sf) + ^ J ik ( r fcifM - r ifci (cj)) . 



(7) 
(8) 



The three operator Green's functions 

r lfc >)=«sfs- ; s'+» w 

must also obey their corresponding equations 

uT ikj (u) = (<% - S kj ) (StS£) u 

+ 1/2 Ja « (S? Si - S^S+) S^;S+ >X 



+ J ki « S? (S?S k - S* k Sf) ; S+ », (9) 
i 

As higher-order Green's functions are generated, one ob- 
tains an infinite hierarchy of equations which have to be 
decoupled in order to obtain Xij (w) . 



In a frustrated lattice, one is constrained to the non- 
magnetic phase, where (5?) = (Sf) — and (S^S^) = 
2<SfSj). 

Here we adopt an extension of the Kondo-Yamaji de- 
coupling approximation for the kagome lattice, assuming 
(only for different sites i,l,k) 

< Sr S+ ; 5+ > w w a(S^Sf)xkj(u) 

+a(S+S i -) X «H (10) 
« S!SfS k ;St « a (SfSf) * fe ». (H) 

The static correlation functions (5*^5") between spins 
at the lattice positions R^ and Rj depend on the dis- 



tance d 



R,- — R, 



on the number n of bonds be- 



tween sites i and j, and on the lattice symmetry. As a 
shorthand notation, we shall call the local correlation as 
Co = (S^S^), the nearest-neighbor correlation as c\ and 
the next-nearest-neighbor correlations as C2 and c' 2 (the 
former corresponding to the shortest distance, as illus- 
trated in figure pi). We also introduce c„ = ac„\ T r 



parameter a is useful to enforce the condition Co = l/2l 

These correlation functions can be evaluated from the 
spin susceptibility as 



i 

7T 



Imxij (w + 477) 



1 



(12) 



Thus, c i, C2, a nd c' 2 and also the parameter a appearing 
in Eqs. |K)| - |ll| must be evaluated selfconsistently. 

Applying approximation |i~o[]ll| to the kagome lattice 
yields 

\u 2 — 4J 2 c] Xij(w) = — 8Jci6ij + 2c\Jij 
+ Jc!^2 J'ikXkj{u)- J (c + 2ci)'J2Ji k Xkj(u), (13) 



where 



c = 1/2 + ci + c 2 + 5 2 



(14) 



and 



j., I J if i, j are next-nearest neighbors, 
1J I otherwise. 



By working directly with the c$ , the numerical prob- 
lem is reduced to only two selfconsistent parameters, c.\ 
and C2 + c' 2 , from which one can obtain a = 2cq. 

With reference to the underlying triangular lattice, the 
approximate equation O can be rewritten in matrix form 



-2ci (iJSijl 



Ji 



2 - 4J 2 c] 

JciY J'ikXkj(u) ~ J (c + 2ci)^ J ik x kj (w) , 

k k 

(16) 
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where J y - and J\ • are the exchange matrices connecting 
neighboring triangles. As indicated in Fig. 1, each tri- 
angle is surrounded by six neighboring triangles, all of 
them contributing to both types of exchange matrices. 

The elements of the dynamic susceptibility matrix are 
defined as 

X° 6 fa> ^ = W ? e_lq ' (R ^ R?) ( 17 ) 

where R" = R^ + T , R^ is the lattice position of the 
i-th triangle, and the vectors T a indicate the position of 
site a in the triangular unit cell (a = 1, 2, 3). 

Applying a corresponding Fourier transformation to 



equation [l6| gives 

X (q,w) = -2c 1 A- 1 (q,^) B(q) (18) 

where 

A(q,w) = (lo 2 - 4 J 2 c) 1 - JciJ'(q) 

+J(5 + 25x) J(q), (19) 
B(q) = 4Ji - J(q), (20) 

where J(q) and «/'(q) are the Fourier transform of the 
exchange matrices Jy and J^, 



J(q) = 2J 



/ 

cos[(x — v / 3y)/4] 
\ cos[(x + \/3y)/4] 



cos[(:r - \/%)/4] cos[(x + V3?/)/4] \ 
cos(x/2) 
cos(a;/2) / 



(21) 



J'(q) = 2J 



/ 2cos(x/2)cos(v / 32//2) cos[(3x + \/3y)/4] cos[(3x - V3j/)/4] \ 

cos[(3x + V / 3y)/4] cosx + cos[(a;-\/3y)/2] cos(-\/32//2) 



V cos[(3x - v / 3y)/4] 



s(v%/2) 



(22) 



cosx + cos 



i[(x + v / 3y)/2] / 



in which a; = q x a and ?/ = q y a. 
It follows that 

X(q,w) = -2ci 



P(qV) 
A(q,^) 



(23) 



where A(q, £) is a third degree polynomial in £ with q- 
dependent coefficients and the elements of P(q, £) are 
second degree polynomials in £ with q-dependent coefi- 
cients. The dynamic susceptibility can be finally writen 
as 



x (q )W )=-2 Cl ^— — 



Hi(q) 



/=i 



fif(q) 



where 



and 



A'(q,fif) 



A'(q,£) = -A(q,£) 



(24) 



(25) 



(26) 



The f2;(q) are excitation energies, obtained from the 
roots & = f2?(q) of the polynomial A(q, £). One of these 



r 



roots, f2f(q) = 6J 2 (c + c~i), is found to be non-dispersive 
over the whole Brillouin zone. Thus, the dispersion in 
the energy spectrum is due to the other roots (/ = 2,3), 
which are evaluated numerically. 

By writing the correlation functions of eq. |l2| in matrix 
form and performing the inverse Fourier transformation, 
we can show that 

3 

4 b = -2d E Mq) Rf(l) cos{q • (R? - Rj)} 

1 = 1 q 

(27) 



where 



2T^ : 



(28) 



and n(w) is the Bose distribution function. The sum over 
momenta must be performed numerically. 

The dynamic structure factor matrix is related to the 
dynamic susceptibility by 



JS(q,w)=- 



21m x(q,w) 



(29) 



1 - e-P" 

After obtaining the selfconsistent parameters, the static 



0.15 




3 

T 



T=0 

T= 1.466 J 
1=1.151 J 



FIG. 3: Correlation functions versus distance (in units of the 
lattice parameter a). 



FIG. 2: (a) The selfconsistent parameter a and (b) the cor- 
relation functions ci, C2 and c' 2 as a function of temperature. 

structure factor can be evaluated by 

3 

S(q) = 2 Cl ]T^Mq) fl,(q). (30) 
;=i q 

III. RESULTS 

We have calculated the selfconsistent values of ci, C2, 
and as a function of temperature. They are plotted in 
fig. |^-b. The selfconsistent parameter a is shown in fig. 

a. From the knowledge of these selfconsistent functions 
all missing correlation functions can be evaluated. In fig- 
ure H we show the correlation functions c£ between sites 
separated by at most 6 lattice bonds. The distribution 
of the points is very similar to that obtained in Rcf. |24| 
for the pyrochlore lattice, indicating an exponential de- 
cay with a correlation length of the order of the lattice 
parameter a. 

The internal energy is given by E = —GNJci. The 
ground-state value E = — 0.859NJ is in good agree- 
ment with earlier calculations (see Refs. |^, ^4^5|) . In fig- 
ure [I] we present the specific heat C = dE/dT. This 
curve has a single peak and is very similar to the result 
of high-temperature expansionslij reproducing the high- 
temperature limit C ~ T~ 2 . In the low-temperature 
limit, the time involved in the numerical computation 
increases with decreasing temperature, because the con- 
vergence of the integral in q becomes slower, and in addi- 



tion one needs a higher precision in order to perform the 
derivative involved in the evaluation of the specific heat. 
Nevertheless, a careful analysis of the numerical results 
indicates unambiguously a T 2 dependence of the specific 
heat which extends to very low temperatures. Therefore, 
there will be no second peak in the present model. The 
fact that the specific heat does not vanish exponentially 
indicates the presence of low energy states even though 
the low-temperature peak found in Ref. [l^ is missing. 
Indeed, it should be emphasized that it is not sure that 
this peak will persist in the thermodynamic limited which 
supports the results obtained in the present work. The 
change in entropy J Q T ct ^, - dT' is also shown in Fig. 
The total change in entropy is 0.54 In 2, corresponding 
to a ground-state entropy 0.46 In 2, of the same order as 
found in Ref. |o| 

The eigenvalues A„ (q) of the structure factor are shown 
in Fig. pL The result is remarkably similar to what—is 
found from the high-temperature expansions of Ref.E3. 
Whithin the present approach one can easily follow the 
effect of temperature on the structure factor for all 
wavevectors even at very low T. We observe that its 
highest eigenvalue is degenerate over the whole Brillouin 
zone. For the pyroclore, the structure factor is a 4 x 4 
matrix, which has-pne additional eigenvalue. In a pertur- 
bative expansionE.il the two highest eigenvalues are found 
to lie very close to each other, one of them being com- 
pletely degenerate while the other one has a weakly lifted 
degeneracy, which is crucial to reproduce the main fea- 
tures of neutron scattering experiments. As a further 
work, it would be interesting to extend the present ap- 
proach to the pyrochlore antiferromagnet. 
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FIG. 4: Specific heat as a function of temperature (solid line). 
The dashed line gives the integrated entropy. 




FIG. 5: The eigenvalues A n (q) of the structure factor along 
some symmetry directions in q-space at T = (solid lines) 
and for an intermediate value of T (dashed lines). 
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